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Abstract 

We generalize Abrikosov-Gor'kov solution of the problem of weakly coupled 
superconductor with impurities on the case of a multiband superconductor 
with arbitrary interband order parameter anisotropy, including interband sign 
reversal of the order parameter. The solution is given in terms of the effective 
(renormalized) coupling matrix and describes not only Tc suppression but also 
renormalization of the superconducting gap basically at all temperatures. In 
many limiting cases we find analytical solutions for the critical temperature 
suppression. We illustrate our results by numerical calculations for two-band 
model systems. 
February 1, 2008 

I. INTRODUCTION 

Recent advance in the field of high-temperature superconductivity, in particularly discov- 
ery of strong anisotropy of the order parameter, stimulated new interest to the old problem 
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of the effect of (magnetic and non-magnetic) impurity scattering on superconductivity with 
high anisotropy. In a number of theoretical papers pubhshed within the last few years qual- 
itatively new phenomena were uncoveredS^. Moreover, detailed experimental studies of the 
effect of impurities in high-temperature superconductors are underway (see e.g. Refs. [7|-p!0| 
and references therein). 

A specific, but representative case of anisotropic superconductivity is multiband super- 
conductivity (e.g., Ref. p|,[TTHl5|), where the order parameter is different in different bands. 



Allen showed in 19780 (see also Ref. |l^) that a superconductor with a general anisotropy 
can be treated within the same mathematical formalism as a multiband superconductor, 
if one expands the order parameter, pairing interaction, and impurity scattering in terms 
of the Fermi surface harmonics. In this paper we derive a general formula, analogous to 
the Abrikosov-Gor'kov formula for isotropic superconductor^!, but valid for an arbitrary 
multiband system. According to the Allen's formalism, this result is easily generalizable to 
superconductivity with arbitrary angular anisotropy. We will also show explicit results for 
various limiting cases to illustrate the physics of the interplay between impurity scattering 
and gap function anisotropy. We will illustrate the results on a model system with strong 
interband anisotropy, namely one where superconductivity in one of two bands is induced 
by interband proximity effect. 

II. GENERAL THEORY 

Following the standard way of including the impurity scattering in the BCS theorylll, one 
writes the equations for the renormalized frequency ujn and order parameter A„ (n is the 
Matsubara index), which completely define the superconductive properties of the system: 

^'^an = ^^n + ^^7^(7a/3 + 7a/3) (la) 

A„„ = A„ + ^^(7„,-7:^) (lb) 
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A^ = 2nT KApnlQpn- (Ic) 

The general form of these equations for strong coupling and general anisotropy in terms of the 
Fermi surface harmonics can be found in Ref. |l^ Note that according to Allen's terminology 
we work in the disjoint representation, where Fermi surface harmonics are defined separately 
for each sheet of the Fermi surface, and take into account only the lowest harmonic for each 
sheet. Other notations in Eqs.l have their usual meaning: uOn = {2n + l)7rT, Qan = 
\J'^an + ^on5 7a/3 = UajsNp is the Scattering rate matrix due to nonmagnetic impurities, and 
7a/3 — ^ap^P is the same for magnetic impurities. The coupling matrix A is defined in the 
same way as Allen's matrix Aq,q,/0, Aaa' = Vaa^^Na' ■ Here Na is the partial density of 
states at the Fermi level in the band a. The scattering potential U and the pairing potential 
^pairing symmetric matrices, while 7,7'^, and A are not. We shall also introduce the 
following useful notations: 

A„ = E ^-/3' ^ = E KNJN, N = Y,N^, (2) 

(3 a a 

where is the total density of states, Aq, are partial electron-phonon coupling constants, 
which define the electron mass renormalization in the band a, and A is the total isotropic 
coupling constant, which enters the BCS and Eliashberg equation for isotropic constant gap 
superconductivity. Analogously, we shall introduce partial scattering rates, 

= E ^al3^ = 7a/3 ± 7a/3, 9al3 = ^a^G^ " (3) 

At temperatures close to one can linearize Eqs. 1 with respect to A. To do so, we 
introduce, as usual, the renormalization function Z and the gap function A' : 

Zan = ^anl^ni ^an ~ ^anl'Z'an- (4) 

^an = 1 + Gll2uJn\ ^^^^ + G'+/2a;„) = A„ + E A:,,„r;„,/2a;„, (5) 

a' 

which we can solve for A' : 
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= J2A^,{S.^>+g^^,/2uJn)-\ (6) 

a' 

Now from Eq.^ follows that 

/3n 7 

For weak (7 ^ 27rTc) and for intermediate (7 <^ ujd) scattering the usual trick with sub- 
tracting the clean limit, g = 0, can be applied, and extending summation to infinity (a useful 
matrix formula is (I + A)~ = I — A(l + A) ), one gets 

/37 "7 

where L = ln(27*u;£)/7rTc). 7* ~ 1.78 is the Euler constant. By introducing the eigensystem 
of g^, g^p = J2j Ra^d^R-yi3, we can express X in terms of the difference between the two 
incomplete gamma-functions {4'{x) = I]„>o('^ + a;)"^): 

X^f, = 27rT, ^ ^n\d,/2){u;n + d,l2)-^R,p = ^ R-}^x{d,lA7iT,)R^p, (9) 

7 n 7 

with x{^) = — ^(1/2 + x), which is the standard definition of the matrix function 

X=x(g+/47rT,). 

This result is analogous to the classic one of Abrikosov and Gor 'kovS (AG), but includes 
arbitrary anisotropy. Now solving (||) for L, we find: 

A„ = ^(A-^i + X„^)-iLA^, (10) 
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which means that now Tc is defined by the effective matrix Ag// = (A~^ +X)^^. As it is well 
known in the multiband superconductivity theorylll, in this case is defined by the usual 
BCS equation, Tc = {2'^*ujd/'^) exp(— 1/Amax), where Amax is the maximal eigenvalue of the 
matrix A (in our case, of the matrix A^ff). As can be seen immediately from Eqs.(p|-p!OD 
and the definition of gais, diagonal nonmagnetic scattering rates 7^/3 have dropped out from 
Eq.|TO|. This is the manifestation of the Anderson theorem for a many band case: intraband 
scattering does not influence Tc (in the considered Born limit). As will be discussed below, 
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this argument works only for the intraband non-magnetic scattering, while all other are, in 
principle, pair-breaking. 

Up to the second order in A (assuming that AX is small), 

Aejj = A - AXA. (11) 

If we recall that A forms the eigenvector of A corresponding to its maximal eigenvalue Ae//, 
we can immediately write the lowest-order correction to Ae//: 

^Kff = -Xlff E A^X^^A,/ E Al (12) 

a/3 a 

In the strong scattering case (7 ^ tu/)) this formalism cannot be used. Instead, one should 
use directly. 

III. CRITICAL TEMPERATURE 
A. Weak scattering 

Let us consider explicitly some interesting limiting cases. For weak scattering (7^/3, 7*^ <C 
Tc) one can use Eq.|12|, and expand ^ 0) = tt'^x/2 and write 

When all A's are equal (isotropic case), the standard Abrikosov-Gor'kov result is recovered: 
STcjTc = —T^J^af^O^a/s ~ ^a/3)/^'^co = — {'^ / ^^T^q) J2 Y : ^^at is, nou-maguetic scattering falls 
out. On the other hand, in anisotropic case only the intraband non-magnetic scattering falls 
out of Eq. ^ as for instance in a two-band case: 

^ 721 - 721 27^2 + 721 + 721 y 

_ 7r[A;(27i\ + Y12 + 712) + AiA2(7i^2 + 72i - 7i2 - 721) + ^1(2112 + I21 + 72i)] 

8Te(A2 + Ai) 

The main point of the AG theoryEl is that 7'' enters equations for to and A with opposite 
signs. That is why the magnetic impurities appear to be pair-breakers, and the non-magnetic 



6TJT^ = -—A 



ones not. The above solution shows exphcitly that in the multiband case of Eqs.l only 
intraband non-magnetic scattering does not influence Tc {'~fai3 drop out). In an interesting 
limit of two bands, in which one band is superconducting and another is not, An 7^ 0, A12 = 
A21 = A22 = it follows from Eq. |1^ that 

STjn = -^(27n + 7i2 + 712), (15) 

where the first term is the usual AG Tc-suppression, and the last two show that the pair- 
breaking influence of the non-superconducting band is the same both for magnetic and 
nonmagnetic scattering. However, the sign of the order parameter, induced in the second 
band, is different: the same for nonmagnetic and the opposite for magnetic scattering (cf. 
the A2 = curves in Fig.^). Such sign reversal is discussed in more detail later in the paper. 

In the next order in Xp^a^ti the additional correction to 6Tc/Tc is {6Tc/Tc)i = — gfr[(7|i — 
72i)Ai2 + (7i2 — 7i2)A2i]/(Aii — A22) (this corresponds to the so-called interband tunneling, 
specific cases of which are considered in the literature^l) . In the limit of 7^?^- = the above 
expression coincides with that derived in Ref. |TB|. Since 712/721 = A12/A21 = N2/N1, the last 



expression can also be written as 

(STJT,), = -^(7?2 - 7i2)A2i/(An - A22). (16) 

Note that if A21 = suppression of Tc is independent of A22, as long as An > A22- It is clearly 
seen, for instance, in the left-hand part of Fig.|I|, where the suppression rate for A21 = and 
various A22 is shown, and is practically independent on A22- For producing this figure we 
have solved Eqs. 1 numerically for two bands, assuming A21 = A12 = 0, 7^^ = 711 = 722 = 0, 
and 7i2 = 721- In full agreement with Eqs. p!5|JT6| is first suppressed linearly with the rate 
7r7i2/8Tco, then, at 712 ~ it starts to deviate from linearity, and as it will be proved later 
in the paper, saturates at some value depending on A22- 

Another important limiting case, also often considered in the literature, is the limit of 
the weak anisotropy. Let us assume that = A + 6Aa, where \SAa\ ^ A. The pair- 
breaking effect of magnetic impurities is then given by the isotropic AG theory, so it is 
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sufficient to consider only non-magnetic scattering. Let us also take, for simplicity, an 
isotropic scattering, g'^^ = 'y{Sai3 — 1). Then Eq.|T3| gives: 



where '-ftot is the total non-magnetic scattering, summed over all bands (or Fermi harmonics). 
Thus in case of weak anisotropy the suppression is given by the AG formula with an 



effective scattering rate 'y^ff = 7* + (5A2/A^)7. This result has often been obtained for the 



eff 

angular gap anisotropy^SJi. 



B. Interband sign reversal of the order parameter 

Returning to the two-band case, we observe that Eq.|10| is invariant with respect to si- 
multaneous change of signs of A12 and A21 and interchange of nondiagonal magnetic and 
nonmagnetic scattering 712 — > 712 ? 721 72i- This remarkable property is the consequence 
of the symmetry of the matrix Ag// with respect to the above transformations with simul- 
taneous reversal of the relative signs of the order parameters Ai, A2. One manifestation of 
this phenomenon is discussed above for induced superconductivity. Another illustration is 
given by the symmetric case An = A22 = Ay and A12 = A21 = A_l. Then it follows from Eq. 
|10| that with sign change of A^ the role of magnetic and nonmagnetic interband scattering is 
completely reversed. Namely, for positive A_l (Ai and A2 have the same signs) only magnetic 
interband scattering suppresses according to (Tco — Tc)/Tc ~ 7r(7^2 + 72i)/8^c- In the op- 
posite case of negative A^ ( Ai and A2 have different signs) the magnetic impurities do not in- 
fluence Tc but the nonmagnetic ones suppress it according to {Tcq—Tc)/Tc ~ vr(7i2+72i)/8Tc. 
The case of arbitrary A22/A11 is shown in Fig.|^, where we show the numerical solution of 
Eqs.l for the same model as we used in Fig.|^: An = 0.5 is fixed, and A22 changes from 
to 0.4. Both magnetic (7^2 7^ 0) 7i2 = 0) and non-magnetic (7^2 = 7i2 7^ 0) impurities 
are considered. In the ffist band the order parameter is suppressed equally by magnetic 
and nonmagnetic impurities: the solid curves in Fig.^ are the same for both cases. The 
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order parameter in the second band has the same absolute value for pure magnetic or for 
pure nonmagnetic scattering, but its sign is different in the two cases. Moreover, even if 
Ai2, A21 7^ 0, but A12, A21 ^ All, A22, there still is a possibility of the interband sign reversal 
of the gap due to magnetic impurities. This happens when nondiagonal elements in the 
effective A matrix in Eq.[Ty become negative: Ai2'^ = A12 + 7rAiiA22(7i2 — 7i2)/8rc0) which 
does happen if 7^2 is sufficiently large. Then the order parameters in different bands have 
different signs, i.e. solution with sign(A^) = — sign(AQ,) corresponds to a minimum energy. 
This sign reversal leads to an interesting effect: if one starts from a pure superconductor 
with weak interband coupling, and suppresses Tc by adding interband magnetic impurity 
scattering, at some critical scattering strength suppression rate drops drastically. The final 
comment to Fig.^] is that it shows either solely magnetic, or solely nonmagnetic scattering. 
When both kinds of scattering are present, the order parameter in the second band is much 
smaller than in the either pure case, and becomes zero when magnetic and nonmagnetic 
scattering are equally strong. Numerical illustration of this effect can be found in Ref. ^ 

This situation is closely analogous to the known case of d-pairing, where isotropic non- 
magnetic impurity scattering leads to an AG suppression, but with a factor of two smaller 
coefficient (cf. Refs. E^ipSf). If we label those parts of the Fermi surface that have positive 



order parameter as 1 and those which have negative order parameter as 2, then only in 
the "interband channel" the nonmagnetic impurities are pair-breaking, while the magnetic 
impurities are pair-breaking only in the "intraband channel" . Correspondingly, the effective 
pair-breaking scattering rate will be 7^;^ +722 + 712 + 721 = iltot + ltot)/'^- Note that isotropic 
magnetic scattering results in exactly as much pair-breaking in terms of as isotropic 
non-magnetic scattering, contrary to the popular misconception that only nonmagnetic im- 
purities are suppressing in rf— wave superconductors. Interestingly, if isotropic magnetic 
and non- magnetic scattering are both present, and have equal strength, suppression rate 
is the same for the s— and d— wave superconductors. If only magnetic scattering is present, 
Tc is suppressed twice faster in an s-superconductor. In fact, most of these statements are 
not specific for the cZ-pairing, but are true for any superconductor with zero average order 
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parameter and non-zero average square for the order parameter. Let us, for example, prove 
that in such a superconductor isotropic magnetic and non-magnetic scatterings both have 
the same effect on Tc. According to Eq. [1^, Tc suppression rate is proportional to 

(Ak^+^,Ak)k,k' = (A^G+)k - (Akr-j^,AkOkk' = (A^ (r+^,)k,)k - (Akr^^,Ak')kk', 

where we used k and k' for indices to emphasize that the formalism is valid both for interband 
or for angular anisotropy. For isotropic scattering, 7k^k' = 7) 7k k' ~ 7**; ^^^^ equation reduces 
to 

(A^)(7 + 7^) - (AkAk')(7 - 7^) = ((A^) + (AkAkO)7^ + ((A^) - (AkAk'))7. 

For isotropic s-wave superconductors, (Ak) = (AkAk') = A^, and the Tc suppression rate 
does not depend on 7. For a superconductor where (AkAk') = 0, (Ak) 7^ 0, a specific case 
of which is a rf-wave superconductor, the suppression rate is proportional to (7 + 7^*), as we 
have conjectured before. 

C. Strong scattering 

Let us now go beyond the weak scattering limit, so that we cannot any more use the 
expansion in XA in Eq.|10|. In accordance with the AG result, the critical temperature 
vanishes at some finite rate of intraband magnetic scattering 7*^ ~ Tcq. The situation 
is qualitatively different with respect to interband scattering. We will show that in the 
strongly anisotropic case of An, A22 ^ ^aj^is the critical temperature does not vanish even in 
the regime of very strong interband scattering. Let us first consider intermediate scattering 
regime vrTc <^ •ya/s ^^d- In this case one still can use Eqs.(|8|-p!QD. Using expansion 

x{x — ^ 00) = log(47*a; + const) 

we obtain that 

A^ff = A - A • log(7*g+/7rr + const) ■ A, 
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which has a particularly simple form for the case we are interested in, 711 



722 



7i2 -712 
-721 721 

Ae// 



A - (712 + 721)"^ log[47*(7i2 + 72i)/vrT,]A ■ g+ ■ A. 



T 



712/(712+721) 



.27* (712 +721). 

where we assumed, to be specific, that An > A22- Solving for Tc, we get 



(18) 



cO 



7lT, 



cO 



■ 712/721 



27* (712 + 721) 



cO 



27* (712 + 721; 



N2/N1 



(19) 



where Ni^2 are the densities of states in the two bands (the last equality appears because 

Eq. 19 gives the limit Tc ^ when '~fai3 00. However, Eq.lO becomes invalid in this 
regime, namely when the interband scattering rate jaiB exceeds the characteristic electronic 
energy scale ujd which is relevant for the Cooper pairing. In this case, we have to go back 
to Eq. 1^. This equation can be solved analytically in an important regime of the isotropic 
superstrong interband scattering, 7q,^ = 7(1 — 5ai3)Np. In this regime, g'^^ = '^{Sa/sN — N^), 
where N = J2a ^a- To handle Eq. ^ we first need to transform the matrix 

{2iuJ^^ + g+p)-' = [(2^„ + ^N)5^f, - ^N^]-' = (2a;„ + 7A^)"'[5a/3 - 7A^/3/(2wn + 7^^)]"' 

to a more tractable form. Expanding the square bracket in series in 7A^^/(2u;„ + 7A^) and 
collecting the appropriate terms, we observe that 

which in the sought limit 7 — > 00 is simply No,/2Nun- Thus 



(20) 



which has the solution 
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A„ = A„A^^ (21) 

n 

where Aq, = X]^o/3 is the mass renormahzation parameter, and the average gap A = 
J2a^a^a/N Satisfies the regular BCS equation with the isotropically averaged couphng: 

A = X-AJ2^, \ = Y: N^XJN. (22) 

So, in the superstrong couphng regime saturates at a hmiting value, which is actually the 
critical temperature calculated in fully isotropic BCS theory. This regime corresponds to the 
so-called Cooper limit investigated previously for the proximity-effect coupled systemslll@. 
Note that the order parameters in the individual bands are nevertheless different, specifically, 
Aq, = AAq/A. This does not mean that the observable zero-temperature gaps are going to 
be different. In fact, they are the same, as discussed in the next Section and illustrated on 
Fig. |. 

We illustrate the above discussion of the Tc suppression by a numerical solution of the 
Eqs.0 for a two-band case with strong interband anisotropy and A12, A21 ^ An, A22, tIi = 
7^2! corresponding to equal densities of states in the two bands Ni = N2. We have chosen 
Ai = 0.2, A12, A21 = and = 712 = 0. The results of calculations of Tc vs 712 are shown 
in Fig.|l| for various values of A2. In accordance with the above analytical results, Tc first 
drops steeply as 712 increases and then saturates at some finite value when 712 ~ cud. The 
saturation value depends on A2 in accordance with Eq.^ The suppression of Tc remains the 
same when nonmagnetic scattering is zero, 712 = 0, but 7J2 is finite, except that the order 
parameter A2,n=o has the sign opposite to that of Ai_„=05 as discussed above and illustrated 
on Fig.|^. Both kinds of impurities, 7^2 and 712, suppress the critical temperature in this 
case according to Eqs.|I^-|22. 



IV. DENSITY OF STATES AND SUPERCONDUCTING GAP 



The discussion of the critical temperature suppression was based on the solutions of the 
linearized equations. To obtain the density of states, the nonlinear Eqs.l should be solved. 
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In the presence of impurities there is no distinct gap, in the sense that the minimal excitation 
energy does not coincide with the maximum in the density of states. The latter is defined 
in terms of A^^,. Namely, the superconducting density of states N{uj) in units of the normal 
density of states at the Fermi level A^o is 



where A'(c<j) is the analytical continuation of A'„. An analytical solution for A^^ is not 
straightforward to obtain; however, some properties of the numerically obtained solutions 
for A^Q, are already illustrated above in Fig.^. Moreover, there are some rigorous statements 
that can be made about A^^,. 

Let us consider again the limit of the superstrong isotropic non-magnetic scattering. We 
have shown above that Tc in this case is reduced to Tc of the equivalent BCS system with 
the isotropic coupling constant. The same statement appears to be true for the gap in the 
excitation spectrum. Indeed, following AG, we can define Uan = ^n/ ^'m^ ^^id, after the usual 
transformation Un —iu, Uan —iua, ^'an ~^ A^(co'), we can write down the multiband 
analog of the Eq.42' of AG: 

A„-"" A.^ ^/r^ • ^'^^ 



In the absence of magnetic impurities we let F^^ = F~^ = 7, and tend 7 ^ 00. Evidently, 
a solution of Eq. in this limit exists only if Ua = up, and correspondingly A^ do not 
depend on a. We conclude that in this limit the reduced density of states is the same in all 
bands, and in fact coincides with that of the isotropic BCS model with the gap determined 



from the nonlinear analog of Eq.^ A = A A J2n SvrT/ ycj^ + A^ , with A = J2a Na^a/N. 

The evolution of the densities of states in a multiband non-magnetic scattering case is 
shown on Fig.|^. Here, we show the results of numerical solution of Eqs. 1 in the weak 
coupling regime with Ai = 0.5, A12, A21, A22 = and 7^^ = 0. Only nonmagnetic interband 
scattering 712 = 721 is included. In the clean limit, the two bands show two different 
excitation gaps. In accordance with earlier calculationsBH, any weak, but finite impurity 
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scattering mixes the pairs in the two bands, so that the first band (with the larger gap, i.e. 
more superconducting) develops a tail in the density of states which extends all the way 
down to the second-band gap. Except for this tail, which consists of the normal excitations 
of the second band, scattered into the first band by impurities, the density of states still 
looks similar to the clean-limit one. Upon the increase of the scattering rate, the low-energy 
tail in the first band density of states grows, and the minimal gap, the gap in the second 
band grows as well. This reflects the fact that larger number of pairs is scattered into the 
second band and induced the interband-proximity-effect superconductivity there. Thus the 
decrease in the critical temperature of the system is accompanied by the increase of the 
minimal gap in the excitation spectrum. 

Next, let us include interband magnetic scattering into Eq.^. Then in the considered 
regime 7 ^ oo we have F^^^ = 7+7*, F"^ = 7 — 7^* and = up = u (we assumed l^^p = Y)- 
As a result, the densities of states in each band are given by N{uj) = Reu/Vu"^ — 1, where 
M is a solution of the equation 



An energy gap corresponds to maximum real solution for u in the interval u < 1 and the 
pair-breaking rate is given by 27**. Thus, with increase of non-magnetic scattering we have 
a crossover from the state with different signs of order parameters in different bands (for 
zero 7) to the isotropic state (for 7 —>■ 00). This isotropic state may be normal, gapless or 
gapped, depending on the value of 27*. Following the AG analysis, we obtain that an energy 
gap at 7 — s> 00 will exist if 7* < exp(— 7r/4)Ao/2, where Aq is the BCS gap at T = 0. This 
case is particularly interesting: since in the isotropic (7 00) limit there is a finite gap 
and a finite positive order parameter in both bands, and in the opposite limit of small 7 the 
order parameter in one band is negative, it is clear that at intermediate values of interband 
scattering 7 a gapless state should be crossed over. The last statement is in agreement with 
the result of Ref. |I| that for an order parameter with a sign change and a nonzero Fermi- 
surface average, a gapless state develops with an increase of impurity concentration, but the 




(24) 
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gap is restored at large concentration of impurities. To illustrate such crossover, induced by 
magnetic scattering, we show in Fig.^ the results of our numerical calculations for a weak 
coupling 2-band model (Ai = 0.5, A12, A21, X22 = and 7^2 = I21 = Tc/2. The total density 
of states is shown at different values of interband nonmagnetic scattering rate 7. According 
to the discussion in section [111 B| , for small 7 both order parameters have different signs. In 
this case 7^* has no pair-breaking effect, and small gap (negative order parameter) is induced 
in the second band by the magnetic scattering. On the contrary, the interband 7 is in this 
situation pair-breaking. The shape of the density of states shows two characteristic peaks, 
one, at about 0.1 Aq, due to this induced energy gap in the second band, and another just 
below Aq from the gap in the first band. 

With the increase of the nonmagnetic scattering rate the order parameter in the second 
band becomes smaller in the absolute value, still remaining negative. The lower peak in the 
density of states gets washed out and the minimal energy gap becomes smaller. When 7 
approaches 7* this small gap vanishes, although there is still a distinguishable peak in the 
density of states coming from the gap in the first band. At larger 7 >> 7* both gaps have 
again the same sign and now it is 7'^, which is pair-breaking. As one again can see from 
Fig.^, a small gap is restored for the last two curves, corresponding to 7 = 2Y = T^, and 
7 = 2O7* = lOTc. Note that at 7 >> Tc, the gap cannot any more be ascribable to any 
of the two bands, but corresponds to a fully isotropic superconductivity, as described by 
Eqs.|6l-|2. 



V. CONCLUSIONS 

In conclusion, we generalize the Abrikosov-Gor'kov solution on the case of arbitrary 
interband anisotropy of the pairing interaction, and arbitrary strength and anisotropy of 
magnetic and/or nonmagnetic impurity scattering. The results are illustrated on model 
two-band systems with interband anisotropy and various kinds of impurity scattering. In 
case of weak scattering, we found an analytic solution, analogous to the isotropic solution 
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of Abrikosov-Gor'kov. For weak anisotropy, this solution yields the critical temperature 
suppression proportional to the mean square variation of the order parameter, the fact 
earlier pointed out by several authors in various special cases. We also prove analytically 
that the superconductivity suppression by isotropic magnetic and isotropic non-magnetic 
impurities is exactly the same when the average order parameter is zero {e.g., in case of 
(i-pairing). We also give an analytical solution for Tc in the two-band model in case of 
intermediate-strength scattering. In case of superstrong scattering we find a solution for 
Tc for arbitrary anisotropy. We also discuss the evolution of the density of states with the 
increase of the impurities concentration (or scattering strength) . 
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S.A. Wolf, W.E. Pickett, S.V. Shulga and S. Tolpygo. This research was in part supported 
by the INTAS grant N 94-2452 (AAG) and by the Office of Naval Research (JIM). 
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FIGURES 



FIG. 1. Tc suppression by interband scattering. Dots on the right axis show the asymptotic 



value of Tc at 7 ^ oo, according to Eq. 22 



FIG. 2. Dependence of the order parameters Ai^„=o and A2,n=o in two bands in a two-band 
model with A12 = A21 = and different ratios A22/A11, on the interband impurity scattering. Solid 
lines show the order parameter in the first ("superconducting"), dashed and dotted lines show the 
order parameter in the second band, where superconductivity is induced by impurity scattering. 
Dashed lines correspond to non-magnetic interband scattering, dotted lines to magnetic interband 
scattering. 

FIG. 3. Superconducting density of states in a two-band model with A12 = A21 = A22 = 0, 
All = 0.5. Only interband non-magnetic scattering is included. The solid lines show DOS in the 
band 1 and the dashed lines in the band 2. Note that both DOS coincide in the regime of strong 
scattering. 

FIG. 4. Total superconducting density of states in a two-band model with A12 = A21 = A22 = 0, 
All = 0.5. 
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